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Abstract. We show that in general, the specification of a contact angle condition at the 
contact line in inviscid fluid motions is incompatible with the classical field equations 
and boundary conditions generally applicable to them. The limited conditions under 
which such a specification is permissible are derived; however, these include cases 
where the static meniscus is not flat. In view of this situation, the status of the many 
'solutions' in the literature which prescribe a contact angle in potential flows comes into 
question. We suggest that these solutions which attempt to incorporate a phenomeno- 
logical, but incompatible, condition are in some, imprecise sense 'weak-type solutions'; 
they satisfy or are likely to satisfy, at least in the limit, the governing equations and 
boundary conditions everywhere except in the neighbourhood of the contact line. We 
discuss the implications of the result for the analysis of inviscid flows with free surfaces. 
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1. Introduction 

Consider an inviscid liquid, under a passive, inert gas, partially filling a smooth walled 
container. A gravitational field acts on the liquid. The liquid-gas interface is subject to 
surface tension. The interface meets the walls of the container at a line called the contact 
line. At any point on the contact line, the angle between the normal to the gas-liquid 
interface and the normal to the solid wall is called the contact angle, a. In the quiescent 
state, for many pure materials and smooth solid surfaces, the contact angle for many 
gas/liquid/solid systems is a function of the materials alone. Here, we will assume this to 
be true and call this contact angle, the static contact angle a s . 

The surface tension at the gas-liquid interface, from now on called the interface, 
requires that there be a jump in the normal stress across it if it is not flat. In the inviscid 
case the jump is in the pressure and for a static interface the liquid pressure is just the 
hydrostatic pressure. Thus the shape of the static interface depends on a Bond number 
Bo, the ratio of a measure of the gravitational force to a measure of the force due to sur- 
face tension. But this shape also depends on the static contact angle a s , which provides a 
boundary condition for the differential equation that determines the static interface shape. 
The static contact angle is therefore a parameter that influences the static meniscus shape 
and needs to be prescribed in order to calculate the meniscus shape. 

When the interface is in motion, the surface tension again requires that a pressure 
jump, proportional to the interface curvature and to the coefficient of surface tension, exist 
across it. Of course, the pressure in the liquid will now be determined by the unsteady 
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Bernoulli equation. The situation as regards the contact angle, however, is more compli- 
cated and confused. If one examines the dynamical equations for the interface it is not at 
all obvious that one needs to prescribe a condition on the contact angle. Moreover, as is 
well-known, for the case of linearized disturbances between plane, vertical walls where 
the initial interface is plane, the classical solution can be obtained without any specifica- 
tion of the contact angle, which turns out to be constant and is equal to n/2 throughout 
the motion. On the other hand, there are many examples in the literature where analyses 
and calculations have been made of unsteady potential motions where a condition on the 
contact angle has been prescribed as a boundary condition at the contact line. Just a short 
list of these could include Miles [9 1, Billingham 1 3 1 and Shankar [ 15 1. It appears that the 
motivation to prescribe the contact angle comes from the apparent behaviour of real, vis- 
cous interfaces and the need to tailor inviscid models so that they lead to realistic results 
for real interfaces. The question that we raise here is: are we really free to prescribe the 
contact angle in inviscid potential flows and if not, what is the status of the 'solutions' 
that have appeared in the literature that purport to model 'real' contact angle behaviour? 

Our interest in this question is recent and followed an investigation of contact angle 
behaviour in viscous flows with pinned contact lines [ 17 1. We had been aware that Ben- 
jamin and Ursell [2 1 had shown that the contact angle would remain constant in linearized, 
inviscid, potential motions about a flat interface, i.e. one corresponding to a contact angle 
of n/2. Without carrying out an analysis, we had assumed the result to be generally true 
and in 1 17 1 had, perhaps influenced by all the work employing a constant contact angle, 
even asserted this. When, however, we recently tried to prove the result, it began to be 
clear that there was a problem here: analysis, following [2 1, seemed to show that the con- 
tact angle cannot, in general, be prescribed. Our purpose here is to demonstrate this and 
to try to place the existing literature (including ours!), in which a contact angle condition 
is employed in potential motions, in proper perspective. We believe that this is an impor- 
tant issue because even if an inviscid 'solution' is a good model of reality in some sense, 
we should be clear in what sense, or approximate sense, the 'solution' is a solution. 

2. Analysis 

We consider the inviscid motion of a liquid in an arbitrary, three-dimensional smooth 
walled container. The motion is generated by the translational motion of the container. 
The restriction to translational motions is essential to ensure potential flow in the moving 
frame, a necessary condition for some of the results that will be derived. The fluid is 
initially in static equilibrium with the gas above it which is at uniform pressure 1 . The 
motion is assumed to start and continue with a uniform pressure over the interface; we 
will assume the gas to be passive, i.e. it only exerts a constant pressure on the liquid 
interface. In §2.1, we write down the equations governing the motion. In §2.2, we first 
consider planar motions; motions in a cylinder of arbitrary cross-section are examined in 
§2.3. Finally, in §2.4, we summarize the main results. 

2.1 Governing equations 

We write the equations in a reference frame attached to the container (please, see fig- 
ure 1). The container wall is given by f(x,y,z) = 0. Rectangular cartesian coordinates 

'The condition of static equilibrium can be relaxed and, in fact, has to be in the case of time periodic wave 
motions. 
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are employed with gravity generally in the negative z-direction. Our analysis is restricted 
to the case when the interface is representable by a single valued smooth function, e.g. 
by z = C(x,y,f). The interface motion can be of finite amplitude however. The equations 
governing the liquid motion are the continuity and Euler equations: 

Vu = 0, (1) 
u, + (u- V)u= -Vp + F, (2) 

where u is the liquid velocity, p is the pressure and F is the net body force, both real 
and fictitious. F can be an arbitrary function of time. These have to be solved subject 
to the boundary conditions on (a) the container wall and (b) the interface. The condition 
on the wall is the no-penetration condition u • V/ = 0. The conditions on the interface 

z=C,{x,y,t) are 

Q + u£x + v£ y = w, (3 a) 

1 

P-/>a=-^-K (3b) 
Bo 

which respectively are the kinematic condition on the interface and the normal stress con- 
dition on it. Bo is a suitably defined Bond number, p a is the constant ambient pressure over 
the interface and K is the local interface curvature. The interface z = £(x,y,t) is assumed 
to intersect the container wall f(x,y,z) = in a smooth curve called the contact line. It 
is further assumed that £, u and p are analytic in all variables and that the smoothness of 
the solutions is upto and including the boundary. It must be pointed out that the existence 
of such solutions is not known at present. The possibly time dependent angle a made by 
the contact line with the wall is given by 

fl-fx^-fy^ =CQsa(f)j (4) 



IV/H1 + C2 + C, 2 ] 1 / 2 

where (4) is to be evaluated at any point (x c (t),y c (t), £(x c (t), y c (t),t)) on the contact 
line. We will show, over the next two sub-sections, that no a priori prescription of 
a(t) is possible though in a few cases, it turns out to be a constant equal to its initial 
value a s . 



230 P N Shankar and R Kidambi 



2.2 Planar motions 

The container is of arbitrary cross-section in the x-z plane and the motion is 2D. The 
contact line in this case consists of just two points (A and B, say). During the course of 
the motion, let the contact line traverse regions R of the wall surface. We will distinguish 
two cases depending on whether R is locally flat or not. We designate the former the 'flat 
or straight-wall' case and the latter as the 'curved-wall' case. Examples of the former are 
rectangular and wedge-shaped containers (with the z-axis along the wedge); a cylindrical 
channel is curved-walled. We first consider the straight-wall case. 

2.2.1 The flat or straight-wall case. The x- and z-axes are chosen such that, locally, the 
body surface is a line of constant x, say x = 0. There will be an x component of gravity in 
this case; the gravity direction has no bearing on the analysis however. The no-penetration 
condition implies w(0,z,f) = from which it follows that k („) (0,z,f ) = where w („) is 
the nth derivative of u with respect to z. 

Now, differentiate (3a) with respect to x to obtain an equation valid on the interface and 
hence on the contact line 

Cxt(x,t) +u(x,^(x,t),t)Cxx{x,t) + [u z (x,C{x,t),t)Cx(x,t) 

+ u x (x,^(x,t),t)]^ x (x,t)^w x (x,^(x,t),t)+w z (x,^(x,t),t)^ x (x,t). (5a) 

Using the fact that an inviscid flow starting from rest has to be irrotational both in an 
inertial frame and a frame translating with respect to the inertial, we get u z — w x = 0, 
which with the continuity equation u x + w z = allows (5 a) to be written in the form 

Co = u z - (2m x + Cxu z ) Cx - « Ca ■ (5b) 

where the arguments in (5 a) have been dropped to reduce the clutter. When (5b) is applied 
at the contact line we obtain 

£*(0,0 = -211,(0, C (0,0,0^(0,0- (6) 

Note that, if the container were to be rotating as well, u z — w x j^0 and the above analysis 
would not apply. Equation (6) shows that at the contact line cannot be specified arbi- 
trarily; it certainly is not zero in general 2 . This means that the contact angle changes with 
time in a manner that cannot be prescribed beforehand. However, for a = 7c/2, not only 
is (0,0) = but also ^« (0,0) = for all k = 2,3, . ... We will show this by mathe- 
matical induction 3 . Let the induction proposition be 

P(k): Q„o(0,0)=OVm = 0,...,L 

P(0) is true because a s = n/2. Assume P{k) is true; we will show that P(k + 1) is true. 
Rewriting (5b) as = u z — a(u, Q£ x — u£ xx , the kth time derivative of this equation can 
be written as 

C«(*+i) = « a (t) - [aQ,(k) - (uCxx) t (k) ■ 



2 It should be noted that we cannot set u x (0, f (Q,t),t) = and argue that the contact angle is constant; the reason 
is that one cannot specify two conditions at the boundary in a potential flow and w(0,z,() has already been set 
toO. 

3 It is not possible to integrate (6) directly in time to establish the result since all that we know is that ^(0, 0) = 
and hence from (6) that £ xt (0, 0) = 0. 
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With respect to the above equation, we observe the following: 

1. d{K z }/df* = 0atx = 0Vifc = 0,l,2,... 

2. d{a£ x }/dt k = Y,m=ob m a t [k-m)£, xt (m) where b m is the binomial coefficient C(k,m). By 
the truth of P(k), this sum is zero at x = 0. 

3. d{u^ xx }/dt k = Ylin=QC{k,m)u t (k- m ) Cxxtim)- This sum is zero Mk = 0, 1, . . . as u and all 
its time derivatives vanish on the contact line. 

Thus we have C xr (*+i) (0,0) = which means P(k+ 1) is true. Thus, by mathematical 
induction, P(n) is true Vn = 0, 1 , . . . and so £ x( (0, t ) = 0. This in turn implies that £ x (0, t ) = 
for all time and the contact angle remains at n/2 for all time. Some observations are 
noteworthy: 

1 . Irrotationality of the motion is necessary but not sufficient. 

2. This is a nonlinear result, i.e. it holds irrespective of the perturbation amplitude and the 
shape of the static meniscus as long as the initial contact angle a s is n/2. In particular, 
the static meniscus need not be flat. 

3. The result holds as long as the region of the body surface over which the contact line 
moves is flat. The shape of the body elsewhere is immaterial. 



2.2.2 The curved wall case. While (5b) is still true, the observations following it are not 
and nothing can be said about the behaviour of the contact angle for nonlinear motions 
even for the case of a s = n/2. It turns out that the contact angle is not constant even for 
linearized motions with a s = n/2. This is counter-intuitive as locally one would expect 
the curved wall to look 'straight' and it is instructive to see where the straight-wall analysis 
breaks down in this case. Representing the body surface by z = g(x) the no-penetration 
condition is written as 

ug — w = 0. (7a) 

Differentiating (7a) on the wall yields 

g"u(x,g{x))+g'{u x + u z g'} =w x + g'w z (7b) 

a s = n/2 yields the relations r/ s '(x t ) = -l/g' and g"(x c ) = r]!//r]^, where T7 s (x) is the 
static meniscus and the primes denote differentiation with respect to x. Using these rela- 
tions in (7b) yields 

7]"u(x c ,z c ) = ri' s 2 w x + ri' s (u x - w z ) - u z . (7c) 

The initial interface may be flat or curved depending on the value of g' (infinite or finite) 
at the contact line. If (3a) is linearized about this initial interface we have 

ri, + uri!, = w, (8a) 

where rj(x,t) is the perturbation, caused by the container motion, of the static meniscus 
77s. We need to show that r] x (x c ,t) = r] x (x c ,0) for all time. To this end, differentiate (8a) 
with respect to x to obtain 



rixt + {u x + ri' s u z )ri' s + ur}" = w x 



(8b) 
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which on using (7c) becomes 

r\xt = (1 - n?)(w x + u z ) + t\' s (w z - 2u x ). (8c) 

Thus TjxtlxoO) = as the fluid is at rest initially. However, the time derivatives of u and 
w are in general not zero which means that rj xt (x c ,t) ^ for arbitrary time t . 

2.3 Three-dimensional motions 

An analysis similar to the one in §2.2 shows that the contact angle cannot be prescribed 
arbitrarily in the three-dimensional case as well. However, following Benjamin and Ursell 
0, we will now show that a contact angle of n/2 remains constant for linear (infinitesimal 
amplitude) wave motions in a right cylinder of arbitrary but smooth cross-section standing 
on one of its ends. Let the body cross-section be given by g(x,y) = and the interface by 
%(x,y,z,t) =z — £(x,y,t) = 0, then we have the body and interface unit normals as 

8x1 + gy] 



whence the contact angle a is given by 

, . . „, dt; 

cos a — n<o ■ n\ = n\, ■ V C = — . 

dn b 

Since a s = 7T/2, we have initially <9£/<9«b = 0. Linearizing (3a) and making use of the 
irrotationality of the motion, the kinematic condition is written as 

where <p is the velocity potential governing the motion. Note that (9) is applied on z = 0. 
Differentiating (9) in the direction of «b, we have 

d 2 <j> 
di%dz 

which on using the no-penetration condition on the contact line d(j>/dn\, = 0, leads to 
C«b' = on the contact line. Since this is true for all time, this means <9£/<9«t, = for all 
time, i.e., the contact angle remains n/2. 

The same result does not hold in the nonlinear case. We show this for the simplest 
case of a right circular cylinder. Employing cylindrical coordinates (r, 8 , z) , the kinematic 
condition is 

£ + M -£ + --i= w (11) 

or r dB 

applied on z = £(r, 0,t) and u,v and w are the r, 9 and z components of velocity. a s = n/2 
translates to d £ / dr = on r = a where a is the radius of the cylinder. Differentiating (11) 
with respect to r, we have 

v v C,q 

Qrt + UQ r r+(u r + U z Q r )£r + -Qre - ~jCe H (>V + V z Cr) = W r + W z £r- (12a) 

r r L r 
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Using irrotationality and continuity, the above equation can be written as 



. (12b) 



Since the motion starts from rest, u — v = w — initially. So are all the spatial derivatives 
of velocity. Initially and £ r g are both zero on the contact line. Finally, £ rt (a, 9,0) = 0. 
For this to hold for all time, we need to show that the higher derivatives vanish as well, as 
we did in the 2D case. Differentiating (12b) once, we have 

Cm = (1 - Cr )(% + u zzQ) ~ 2« z CrCrt ~ 2 ("n + U rz Q)Cr ~ 2-UrCrt 
- {Ut + ll Z Q)Crr - u£rn + ~ (— - CrOl) + ~ (— - Cr9 



r \ r I r \ r 

x ( v t + v zQ) —( v r + V z Cr) ~i v tr + V n Q + ^Crt 

r r 

+ £r(% + v zz Q)} - — v e - —{v e , + v dz Q). (12c) 
r r 

It can be shown that £ rff (a, 0,0) = 0. However the time derivative of (12c) will contain 
terms like ^g t v tr /r which are not necessarily zero and hence £,,(a, 9,t) ^ for arbitrary 
time t . It can be shown however that a contact angle of n/2 is preserved by the class of 
axisymmetric motions (see Appendix C). This is consistent with the result obtained in the 
two-dimensional case. 



2.4 A summary of the results 



We summarize, for the convenience of the readers, the main points of the last two sections. 
The most important point is that in inviscid fluid motions starting from rest in a container, 
the contact angle cannot be prescribed in advance and neither is there need for such a pre- 
scription. However, if one insists on prescribing the contact angle, this would necessarily 
result in a 'weak-type' solution - one that is in violation of the actual behaviour of the 
contact angle. However, in the special case of a s = 7f/2, there exist situations where the 
contact angle remains constant throughout the motion. These situations include cases of 
curved static menisci. Prescription of conditions other than this in this case will again lead 
to 'weak- type' solutions. 



3. Discussion 

Our discussion will center on the considerable confusion that exists in the literature for 
the last five decades, on the role of the contact angle in inviscid fluid motions. First we 
should make clear that there is no confusion whatsoever in the classical literature on wave 
motion in liquids, for example as given in Lamb 1 8 1. In the classical literature, capillarity 
is considered only in situations where the liquids are not bounded by solid boundaries; in 
such situations there is no contact line and so the difficulty we are considering does not 
arise. 
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3.1 Examples from the literature 

We will substantiate our statements above by presenting below a small sample of the 
literature dealing with inviscid contact lines. We hasten to add that we have no wish to 
criticize any particular worker or group; indeed we cite our own work as manisfestations 
of this confused state of affairs. The following is in rough chronological order: 

(a) The article by Reynolds and Satterlee 1 14| appears in NASA SP-106, which was the 
bible for many aerospace engineers for almost three decades. This excellent article 
deals with all aspects of the low-gravity behaviour of liquid propellants. In dealing 
with low-gravity sloshing, they indicate the boundary conditions to be imposed, 'plus 
a contact angle condition' . For the linearized problem they suggest that this takes the 
form 'h r = yK where h is the perturbation to the interface height and 7 is a constant to 
account for 'contact angle hysteresis'. After the formulation of the general linearized 
problem, they indicate that they are still faced with a 'formidable problem'; then they 
just deal with the special flat interface case in a cylindrical tank. For this special case 
they point out "Note that we were unable to enforce any contact point condition" and 
further note that the contact angle is unchanged for this solution. It is to be remarked 
that the authors wished to impose a contact angle condition but were unable to do 
so and found the contact angle to remain constant (= n/2). The reason is that they 
were working with the 'classical' spatial eigenfunctions and could only recover the 
'classical' solution, which our analysis showed will lead to a constant contact angle 
in the flat interface case. 

(b) The papers by Moore and Perko 1101 and Perko | 12 1 are important not only because 
they are among the first to deal with large-amplitude motions and surface instabilities 
leading to breakers, but also because they suggest new methods of dealing with the 
liquid sloshing problem. Both papers deal with the initial value problem for curved 
interfaces under the influence of capillarity and are based on expanding the velocity 
potential in a series of harmonic functions with time varying coefficients. The solu- 
tion method in [ 10 1 is such that the evolution of the interface, at each increment of 
time, does not involve the current contact angle. In fact the only role of the contact 
angle is in determining the initial interface shape. Consequently, we find from their 
interesting figure 3, which shows breakers at the wall, that the contact angle appears 
to change with time. These developments are entirely consistent with the analysis in 
§2 which showed that if a s ^ 7C/2, the contact angle will vary in general. Perko 1 12 1 
extends the earlier analysis to the general axisymmetric case. However, the author also 
points out that it "includes the constant contact angle boundary condition necessary to 
have a well-posed problem in computation" (emphasis added). This is puzzling since 
computations were possible in the earlier paper without any condition on the contact 
angle, and in fact it was not possible to impose such a condition. The author points 
out that the new figure 4, with a constant contact angle, is for the same conditions as 
for the earlier figure 8, where it was not constant; there is no suggestion as to which 
one we should prefer or why. It appears from §2 that the 'solutions' in Perko 1 12 1 are 
not classical ones. 

(c) Chu |4 | is representative of a large class of papers in the 1960s and 1970s which sug- 
gested ways in which slosh frequencies, forces and moments could be calculated for 
axisymmetric containers under low-gravity conditions. This particular paper suggests 
the use of a Galerkin type of procedure. As regards the contact line, Chu says "In 
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addition, there is an interface contact point condition which takes the form ..." and 
gives the condition given in (a) above. 

(d) It was pointed out in §2 that if the initial contact angle is Tt/2 and the side walls are 
straight, even a non-linear, two-dimensional motion would maintain the contact angle 
at n/2. This implies that if non-linear, capillary-gravity standing waves exist, every 
such wave would be a solution of a problem where straight walls are located at the 
nodal points of the wave motion. An example of this possibility is the solution found 
in Concus |5j where such a solution is found to be third order in the amplitude of the 
waves. 

(e) Myshkis et al II II is an encyclopaedic book on low-gravity fluid mechanics. However, 
in Chapter 5 when they formulate the small oscillation problem, after writing down 
eq. (5.2.14) all they have to say on the contact angle is ". . . (5.2.14) is the linearized 
condition for the conservation of the contact angle". In fact, the classical theory tells 
us it will not be conserved in general. 

(f) A totally different aspect is presented by Benjamin and Scott [1|, who appear to 
have been the first to consider the natural frequencies of a confined liquid with a 
flat interface, but whose contact line is pinned. The inviscid modelling of this situ- 
ation is immediately seen to be problematic because the no-slip condition seems to 
be required at the contact line, a condition that cannot in general be satisfied by a 
potential flow. In other words, there is no classical solution to this problem. In light 
of this, Benjamin and Scott using the framework of functional analysis, define vari- 
ous function spaces, operators and other tools to formulate a 'weak solution' to the 
key equation (7b) of their paper and finally get estimates for the frequencies using 
Rayleigh's principle. Moreover, they show that their theoretical estimates agree well 
with measured values of the wave periods. 

(g) Hocking Q is a widely cited paper because a new model for the contact angle 
condition, 'drj'/dt' = X 'dr\' ' jdn\ is introduced for the a s = n/2 case; here r\' 
is the surface elevation of the small disturbance above the flat static interface. 
The model is introduced to account in some way for the 'wetting property' of the 
fluid and includes both the free and pinned edge conditions as limiting cases. In 
his introduction Hocking says "The presence of capillarity adds an extra term to 
the free-surface pressure condition The increase in the order of the pres- 
sure condition, however, requires extra conditions to be imposed when the solution 
is sought in a finite region" (emphasis added). This is a common misconception. 
While it is true that capillarity increases the order of the equation governing the 
static meniscus and hence the number of boundary conditions needed in this case, 
this does not apply in the dynamical situation. This can be seen by just consid- 
ering the flat interface case between vertical walls. In any case, the analysis of 
§2 shows that no extra contact angle condition can be prescribed for classical 
solutions. 

(h) It suffices now to mention that the old confusions persist into the new millenium, 
typical samples being I3I15I16I19I . 



3.2 The current status of the contact angle in inviscid flows and how it has come about 

We will now summarize how the dynamical, inviscid contact angle appears to be viewed 
by most workers and why it has come to be this way. Recall that in the classical works the 



236 P N Shankar and R Kidambi 



question of the contact angle never arose because confined flows with boundaries and cap- 
illarity were not normally studied. Early studies on the latter were confined to linearized, 
two-dimensional flows between vertical walls; here the classical, exact solutions corre- 
spond to a contact angle of 7z/2, which is maintained throughout the motion. In the 1960s 
the space programmes required solutions to more general problems involving curved static 
interfaces and static contact angles other than n/2. The difficulty posed by these problems 
forced approaches that were either semi-analytical or numerical and some like Moore and 
Perko 1 10 1 did not attempt to impose a dynamical contact angle condition; a s affected the 
initial conditions alone through the static meniscus. The imperative to impose a contact 
angle condition at the contact line, not permitted in general by the classical inviscid for- 
mulation, appears to have come from experimental observations of real, viscous contact 
lines. It is well-known I6 116I13I18I that real, dynamic viscous contact lines display com- 
plex behaviour and are not at all well-understood with many parameters playing a role. 
It is in attempting to model this complicated behaviour in an inviscid framework that the 
need for contact angle conditions began to be felt and then applied. The earlier models of 
a constant contact angle and the one used in 1 14 1 (essentially to model contact angle hys- 
teresis) are in a sense non-dynamic. Hocking's |7 1 model is a dynamic one, attempting in 
an inviscid framework to account for contact line hysteresis or viscous wetting effects at 
the contact line. In any case, the purpose is to account for viscous and other real effects in 
an inviscid, potential model of the flow. 

Thus the need to more realistically model the dynamic contact line appears to require 
the freedom to impose a condition at the inviscid contact line. But as was shown in §2 the 
classical field equations and boundary conditions do not in general provide this freedom. 
Then the natural question is, what is the status of the very large and important body of 
work in which a contact angle condition is imposed, in violation of the classical formula- 
tion? Let us call solutions which are obtained without such a condition 'classical'. Then 
this body of work referred to does not deal with classical solutions. This means that these 
'solutions' will be found to violate at least some of the boundary conditions or assump- 
tions at the contact line. 

3.3 The importance of the present result for inviscid free-surface flows 

A natural question that arises is: how important is the present result that the contact angle 
cannot, in general, be prescribed in a potential flow? Before we attempt to answer this 
question we would like to consider the situations shown in figure 2. 

The configurations in (a) and (b) are planar while (c) is axisymmetric and a s = n/2 
in all three cases. According to our theory, finite amplitude motions in (a) preserve the 
contact angle, while the contact angle is not preserved in (b) and (c) even for infinitesimal 
motions. This is not a result that is obvious at all and shows that we must be very careful 
when dealing with the contact angle in inviscid flows. 

Both from general considerations and from the above example it should be clear that 
the result is of real importance in basic inviscid flow theory. But one may still inquire 
whether the result is of any practical importance, for example, in the calculation of the 
natural frequencies of surface waves or in sloshing calculations. This is not easy to answer 
because we do not know how to calculate such things without making any assumptions 
about the contact angle. For example, say we wish to find the natural frequency in a case 
where the static contact angle a ^ n/2. It appears that the frequency cannot come out of 
an eigenvalue problem by assuming a solution harmonic in time and an expansion in terms 
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(a) (b) 




(c) 

Figure 2. The figure shows three examples of liquids in containers which have 

gas-liquid interfaces above them. In all the three cases the static contact angle 
a s = 7c/2, The geometry is planar in (a) and (b), while in (c) the container is a 
cone with the same half angle as in (a), which too is fiat walled. The present theory 
predicts that the contact angle will be preserved in finite amplitude motions in (a) 
but will not be preserved even in infinitesimal motions in (b) and (c). 



of spatial modes 0(x) because the latter would imply some assumption about the contact 
angle at the boundary. It appears that the only way out would be to do the initial value 
problem with the correct static meniscus and let the field evolve freely. But it is unlikely 
that this can be done easily and certainly not analytically. However, we should point out 
that Moore and Perko 1 10 1 solve the initial value problem without making any assumption 
on the contact angle and apparently without making any assumption on a functional form 
for the interface. On the other hand, Perko |12| solves the same problem holding the 
contact angle fixed at the static value. His figure 4 can be compared directly with the 
earlier figure 8 of Moore and Perko. While there are differences, fortunately and as might 
be expected, the qualitative overall pictures are similar. This suggests that the violation of 
the correct condition at the contact line will lead to violation of the interface conditions at 
the contact line and except in extraordinary circumstances will not have much of an effect 
on the overall field. Thus it appears that natural frequencies and sloshing modes will not 
be greatly affected provided the containers are sufficiently large. This must be especially 
true of the natural frequencies which are related to integrals over the whole field. Indirect 
evidence in support of this position is that most frequency calculations done by different 
methods, which presumably violate the correct conditions differently, agree well with one 
another. 
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4. Conclusion 

We have shown that the classical field equations governing the motion of a confined invis- 
cid liquid under a passive gas do not, in general, permit the independent specification of 
a contact angle condition at the contact line. In fact the only cases where such a condition 
may be permissible are when the static contact angle is n/2, the container walls are flat, 
at least in the neighbourhood of the contact line, and (i) the motion is two-dimensional or 
(ii) the motion is a small three-dimensional disturbance from a flat initial interface. The 
restrictions are indeed surprising as is the difference between two- and three-dimensional 
motions. 

These results have a somewhat serious bearing on the vast literature that exists in which 
'solutions' have been found to inviscid motions in which various contact angle conditions 
have been imposed. It is our contention that these cannot be classical solutions to the 
classical field equations since classical solutions do not permit the imposition of a contact 
angle condition. It is suggested that these 'solutions' belong to an improperly defined 
class of 'weak-type solutions', in the sense that they attempt to solve the field equations 
in an approximate sense, with some of the equations being solved exactly. The need for 
such 'solutions' is driven by the compulsion to try to model in an inviscid framework, the 
complicated behaviour of moving viscous contact lines. Examples were given of other 
cases where a similar situation exists. 

Finally, we have shown that while the present result is of basic importance in the theory 
of inviscid free-surface flows, it is unlikely to seriously affect the practical calculation of 
natural frequencies and sloshing modes in containers. 



Appendix A 

In §2, the initial contact angle a s = n/2 was shown to remain constant under certain con- 
ditions. One of these conditions was that the initial interface be the static meniscus itself 
and that the fluid be quiescent initially, i.e., we start from a static equilibrium. This would 
be the normal procedure of posing an initial value problem. On the other hand, we can 
seek special solutions such as periodic motions where though initially the fluid is at rest 
and the pressure over the interface constant, the state is not one of static equilibrium and 
the interface is not the static meniscus. Examples are the well-known linear and nonlinear 
periodic oscillations [8 5| between two parallel vertical walls. The initial contact angle 
a s = n/2 would remain constant in this case as well. 



Appendix B 

We consider here the case of a s = 0; at first glance, it seemed the contact angle would be 
preserved in this case as well. Since r]^(0,t) — °°, we reorient the coordinates such that 
x is along the wall, z normal to it and z = C( x ? f ) — 1 lt( x ) + T ?(*)0 gives the interface, as 
before. The initial interface, close to z = can also be described as x — % s (z), a description 
that will be needed below. 

The linearized kinematic condition is written as 



f]tx = (1 - n?)u z - r\"u - t]' s u x , 



(Bl) 
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where Tjg(O) = 0. 77^(0,0) is zero as the initial state is one of rest. The fcth derivative of 
(Bl)is ' 

rj^+i) = (1 - J? s /2 )M zfW - t] s "m (W - ri' sUxt(k) . (B2) 

It is easy to see that the first and third terms vanish at x = as u z = r]!, — there. If we 
can show that the second term vanishes as well, we would be done. In the second term, 
w(0,0) ^ so we want 77," (0) = 0. 

It is convenient to write the equation for the static meniscus with gravity along the x- 
axis. The equation would be 

1 t" 

e Ss =A, (B3) 



* BO (1 + £*)3/2 

where the prime denotes differentiation with respect to z and A is a non-zero constant 
equal to —k(0)/Bo where fc(0) is the curvature of the interface at the contact line and is 
given by 

K (0) = ^ /(0) /9 - ■ (B4) 

Now, we can write, by chain differentiation, 

*M~f5f§- (B5) 

Since £'(0) is infinite but Jc(0) finite and non-zero, we have from (B4) that £ s "(0)/<i; s /3 (0) 
is a non-zero finite quantity. But by (B5), this means 77" (0) is non-zero as well. Hence, 
^(fc+i) 7^ in general and the contact angle will vary during the motion. 



Appendix C 

In §2, it was shown that the contact angle will, in general, change for nonlinear motions 
in a right circular cylinder. However, we prove here that the contact angle a s = n/2 is 
preserved for the special class of axisymmetric nonlinear motions in a right circular cylin- 
der. The proof by induction follows the general pattern of the proof in §2.2. 1 . Noting that 
axisymmetric motions mean that the interface is given by z = C( r > f ); l et me induction 
proposition be 

P(k): C rtW (0,0)=0Vm = 0,...,jfc. 

P(0) is true because a s = n/2. Assume P(k) is true; we will show that P(k+ 1) is true. 
Writing (12b) for this case, we have 

Q = -ut; rr -2u r t; r +u z (i-tf) (ci) 

as the rest of the terms are identically zero due to axisymmetry. The kth derivative of (CI) 
will have the following terms: 

1. d{u^rr}/dt k — Y,m=oC(k, m )u t (k-m) C rrt W • This sum is zero Mk = 0, 1, . . . as u and all 
its time derivatives vanish on the contact line. 
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2. d{u r £ r }/dt k = L^, =0 C(A:,OT)M rt (j_ m )^ rt ( m ). By the truth of P(k), this sum is zero at 
r = a. 

3. d{u z (l - Q))jdt k = Y l k m=0 C(k,m)u zt{k ^ m) (l - Q) t{m) . This sum is zero VA: = 0,1,... 
as u z and all its time derivatives vanish on the contact line. 

Thus we have C rt (t+i) (0, 0) = which means P(k+ 1) is true. Thus, by mathematical 
induction, P(n) is true Vn = 0, 1, . . . and so £ rt (0,t) = 0. This in turn implies that £ r (0, t ) = 
for all time and the contact angle remains aX.n/2 for all time. 



References 

[1] Benjamin T B and Scott J C, Gravity-capillary waves with edge constraints, /. Fluid 

Mech. 92 (1979) 241-267 
[2] Benjamin T B and Ursell F, The stability of the plane free surface of a liquid in vertical 

periodic motion, Proc. R. Soc. London A225 (1954) 505-515 
[3] Billingham J, Nonlinear sloshing in zero gravity, J. Fluid Mech. 464 (2002) 365-391 
[4] Chu W-H, Low-gravity fuel sloshing in an arbitrary axisymmetric rigid tank, ASME 

J. Appl. Mech. 37 (1970) 828-837 
[5] Concus P, Standing capillary-gravity waves of finite amplitude, /. Fluid Mech. 14 (1962) 

568-576 

[6] Dussan V E B, On the spreading of liquids on solid surfaces: static and dynamic contact 

lines, Annu. Rev. Fluid Mech. 11 (1979) 371-400 
[7] Hocking L M, The damping of capillary-gravity waves at a rigid boundary, /. Fluid 

Mech. 179 (1987) 253-266 
[8] Lamb H, Hydrodynamics (Dover) (1945) 

[9] Miles J, The capillary boundary layer for standing waves, J. Fluid Mech. 222 (1991) 
197-205 

[10] Moore R E and Perko L M, Inviscid fluid flow in an accelerating container, /. Fluid 

Mech. 22 (1965) 305-320 
[11] Myshkis A D, Babskii V G, Kopachevskii N D, Slobozhanin L A and Tyuptsov A D, 

Low-gravity fluid mechanics (Springer- Verlag) (1987) 
[12] Perko L M, Large-amplitude motions of a liquid- vapour interface in an accelerating 

container, /. Fluid Mech. 35 (1969) 77-96 
[13] Pomeau Y, Recent progress in the moving contact line problem: a review, C.R. 

Mecanique 330 (2002) 207-222 
[14] Reynolds W C and Satterlee H M, Liquid propellant behaviour at low and zero g, in: 

The Dynamic Behaviour of Liquids in Moving Containers, (ed.) Abramson H N (NASA 

publication SP-106) (1966) Chapter 11, pp. 387^139 
[15] Shankar P N, A simple method for studying low-gravity sloshing frequencies, Proc. R. 

Soc. A459 (2003) 3109-3130 
[16] Shankar P N and Kidambi R, A modal method for finite amplitude, nonlinear sloshing, 

Pramana - J. Phys. 59 (2002) 63 1-65 1 
[17] Shankar P N and Kidambi R, The constancy of the contact angle in viscous wave motions 

with pinned contact lines, NAL Project Document (2002) CF 0214. 
[18] Shikhmurzaev Y D, Moving contact lines in liquid/liquid/solid systems, /. Fluid Mech. 

334 (1997) 211-249 

[19] Utsumi M, Low-gravity sloshing in an axisymmetrical container excited in the axial 
direction, ASME J. Appl. Mech. 67 (2000) 344-354 



